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O '. Abstract 

We study the role of non-perturbative quantum gravity effects in the Ekpyrotic/Cyclic model 
using the effective framework of loop quantum cosmology in the presence of anisotropies. We 
show that quantum geometric modifications to the dynamical equations near the Planck scale as 
understood in the quantization of Bianchi-I spacetime in loop quantum cosmology lead to the 
resolution of classical singularity and result in a non-singular transition of the universe from the 
contracting to the expanding branch. In the Planck regime, the universe undergoes multiple small 
bounces and the anisotropic shear remains bounded throughout the evolution. A novel feature, 
which is absent for isotropic models, is a natural turn around of the moduli field from the negative 
region of the potential leading to a cyclic phenomena as envisioned in the original paradigm. Our 
work suggests that incorporation of quantum gravitational effects in the Ekpyrotic/Cyclic model 
may lead to a viable scenario without any violation of the null energy condition. 
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I. INTRODUCTION 



One of the most intriguing issues in cosmology concerns with the state of the universe 
in the earliest epoch of its evolution. It is this phase where answers to some of the most 
difficult questions in conventional cosmology are hidden which includes understanding the 
initial conditions of our universe. It is difficult to understand the latter via cosmological 
models based on classical general relativity (GR) since it break downs at very large spacetime 
curvatures and predicts an initial singularity. Therefore, any viable description of the history 
of our universe must necessarily have inputs from a framework beyond GR. Not surprisingly, 
one expects that any such model must capture quantum aspects of gravity in order to be 
successful. Though a treatment based on a full theory of quantum gravity is beyond the 
scope at the present stage, nevertheless useful insights have been obtained on the nature 
of the universe at Planck scale in various frameworks. In one of the earliest works in 
this direction, Wheeler showed that the spacetime near the initial singularity resembles a 
quantum foam [l| (a picture shared with the analysis of quantization of conformal modes 
in gravity (2|). In recent years a different picture seems to emerge in some models based on 
non-perturbative canonical quantum gravity. In loop quantum cosmology (see Ref . 0] for an 
introductory review) a quantization of symmetry reduced spacetimes based on background 
independent loop quantum gravity, backward evolution of states for different stress-energy 
contents, show existence of a bounce near the Planck scale to a classical contracting branch 
[H, 0, 0] . The quantum bounce is non-singular and allows a unitary evolution across the "big 
bang" without any violation of null energy condition (NEC). These results indicate that the 
universe may have a "pre big bang" branch which plays an important role in the history 
and the fate of the universe. 

Interestingly, though loop quantum cosmology predicts existence of a contracting universe 
preceding ours, the idea of a pre big bang branch of the universe is not new. Various 
paradigms in theoretical cosmology have conjectured an existence of such a pre big bang 
branch and existence of a non-singular bounce is the most crucial element for them to 
be viable. As an example, pre big bang branch is envisioned in Ekpyrotic/Cyclic model 
paradigm 0, HI, @] which is also considered as an alternative to the inflationary scenarios (see 
Ref. [l0( for a review) . It is based on the inputs from M-theory where it is envisioned that the 
universe undergoes cycles of expansion and contraction governed by the inter-brane dynamics 
of two boundary branes in a five dimensional bulk. 1 Our observable universe is hypothesized 
to be confined on one of the boundary branes which interacts only gravitationally with 
the other brane. Attraction between the branes leads to their collision, an event which 
corresponds to a big bang or big crunch singularity for the observable universe. 

The transition from the contracting to the expanding branch when the branes collide and 
separate is a tricky issue in the Ekpyrotic/Cyclic models. Thou gh insig hts on this problem 



from the bulk spacetime perspective have been gained (see Refs. 13|, [14], |l5|), its solution re- 
mains elusive in the 4-dimensional cosmological picture. In its absence a viable non-singular 
cosmological model of the early universe is difficult to construct. Insights on obtaining a non- 
singular transition in the Ekpyrotic/Cyclic model is expected from understanding the role of 
non-perturbative quantum gravitational effects which are well understood in loop quantum 
gravity. In particular a lot of progress has been made in recent years on understanding reso- 



1 Based on string theory, idea of a bouncing universe has also been extensively studied in Pre-Big Bang 
models and more recently in string gas cosmologies ■ 
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lution of singularities in the framework of loop quantum cosmology [3,[l6, 17, 18]. This leads 
to a natural question (which is the focus of this work): Is it possible to obtain a non-singular 
transition in the Ekpyrotic/Cyclic model using inputs from loop quantum cosmology? 

Loop quantum cosmology is a canonical non-perturbative quantization of homogeneous 
spacetimes based on Dirac's method of constraint quantization. The gravitational part 
of the classical phase space is labeled by the Ashtekar variables: SU(2) connection A l a 
(which in a cosmological setting is proportional to rate of change of scale factor at the 
classical level) and conjugate triad E? (proportional to the square of the scale factor). The 
elementary variables used for quantization in loop quantum cosmology are the holonomies of 
the connection components and triads. Quantization proceeds with expressing the classical 
Hamiltonian constraint, consisting of the field strength tensor of the Ashtekar connection, 
in terms of elementary variables. The resulting quantum constraint captures the underlying 
discreteness in quantum geometry (one of the predictions of loop quantum gravity) and turns 
out to be uniformly discrete in volume. This has been shown to be a property of isotropic 
models (for all values of curvature index k) 0, 

H H S3, S 13 and Bianchi-I spacetimes 
(24j where one has rigorous analytical control and extensive numerical simulations have been 
performed. 

Some of the novel features of loop quantization can be described as follows. Let us 
consider a isotropic universe filled with a massless scalar field. Take a state which is peaked 
at the classical trajectory at late times and evolve it with the loop quantum constraint 
towards the big bang. It turns out that the state remains peaked at the classical trajectory 
till spacetime curvature (R) is approximately 1% of the Planck value, however on further 
evolution there are significant departures from the classical theory. Instead of following 
the classical trajectory till the big bang as in the Wheeler-DeWitt theory, evolution leads 
to a bounce of the universe when energy density of the universe reaches a critical value, 
Pcrit = 0.41ppi 0. After the bounce the state then evolves further to peak on the classical 
trajectory for a contracting universe. These turn out to be robust features of the theory. 
Using an exactly solvable model it can be proved that the bounce is a property for all the 
states in the physical Hilbert space and there exists a supremum for the energy density 



operator given by the above value [231 ] . Moreover, the fluctuations are bounded across the 
bounce and the state which is semi-classical at late times post bounce is also semi-classical 
at early times pre bounce [25]. It is important to note that existence of bounce does not 
require any fine tuning of initial conditions or a special choice of matter. In fact for the 
above case, matter satisfies stiff equation of state throughout the evolution. Existence of 
bounce and non-singular evolution has also been shown for scalar field with cosmolo gica l 
constant (for positive 26] as well as negative values [27| ) and the inflationary potential [28]. 
Another important feature of loop quantum cosmology is that it turns out that there exists 
a unique consistent loop quantization for the isotropic spacetimes as well as the Bianchi-I 
spacetime jl^, 3(3 . 



Interestingly, it is possible to write an effective Hamiltonian in loop quantum cosmology 
with a resulting dynamics which approximates the underlying quantum dynamics to an 
excellent accuracy. The effective Hamiltonian can be derived using coherent state techniques 



for different matter sources [31|, [32|, [33J. Assuming that effective dynamics is valid for 
arbitrary matter (an assumption which turns out to be true for various cases), it is possible to 
prove that isotropic flat loop quantum cosmology is generically non-singular and geodesically 
complete [13]. It resolves all known types of cosmological singularities. Similar results 
are expected for curved and anisotropic models. Results obtained in the full quantum 
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theory, numerical simulations and the effective theory strongly indicate that resolution of 
cosmological singularities is natural in loop quantum cosmology. 

The availability of an effective description leads us to explore a viable non-singular bounc- 
ing model in the effective description of loop quantum cosmology with the Ekpyrotic/Cyclic 
model potential. This serves as a first step to investigate the role of non-perturbative quan- 
tum gravity effects as derived in loop quantum gravity to the Ekpyrotic/Cyclic models. 
For the case when anisotropies are absent, such an analysis was performed earlier [35| . 
It was shown that a non-singular bounce is generically possible for the flat isotropic Ekpy- 
rotic/Cyclic model in loop quantum cosmology. Despite this a viable model was not possible 
due to lack of a turn around of the moduli field in the epoch when the branes collide and 
separate, inhibiting a cyclic phenomena. It was shown that unless the Ekpyrotic/Cyclic 
model potential becomes positive when the singularity is approached it is not possible for 
the moduli field to turn around in the process of transition from contracting to the expanding 
branch (or vice versa) and lead to cycles. 2 

However in a realistic universe anisotropies are always present, even if their strength be 
very small. Hence it is important to analyze whether the limitations pointed out in the 



analysis of Ref. 35] were artifacts of the assumption of pure isotropy. It is pertinent to ask 
whether there exists a viable non-singular Ekpyrotic/Cyclic model for effective loop quantum 
dynamics incorporating anisotropic properties of spacetime? As we will show the answer 
turns out to be positive. We will see that the mere presence of shear term in the cosmological 
dynamics leads to a turn around of the moduli field from the negative region of the potential. 
Instead of a single bounce of the scale factor, the anisotropic effective dynamics obtained 
from the loop quantization of Bianchi-I model exhibits bounces for each of the directional 
scale factors. Numerical simulations show that the turn around of the moduli occurs in 
the middle of the transition of the mean scale factor from the contracting to the expanding 
branch. Further, the anisotropic shear remains bounded throughout the evolution. Thus a 
non-singular transition from the contracting to the expanding phase with a turn around of 
the moduli as envisioned in the Ekpyrotic/Cyclic model paradigm is achieved in the effective 
loop quantum dynamics. This is the main result of our analysis. 

This paper is organized as follows. In the next section we re-visit the classical theory of 
Bianchi-I model in the Ashtekar variables. We describe the way classical generalized Fried- 
man equation for an anisotropic model can be derived in a Hamiltonian treatment. In Sec. 
Ill we consider the effective Hamiltonian of the Bianchi-I model for the loop quantization 



performed in Ref. 24]. This quantization for Bianchi-I model has recently been shown to 



be the only consistent choice which leads to a bounded shear and expansion factors [30 
Using Hamilton's equations we derive the dynamical equations for connections and triad 
components (the dynamical equations for matter such as the Klein-Gordon equation are 
not changed by quantum gravitational effects) and highlight some of the properties of the 
effective theory. In Sec IV, we use the effective dynamics to analyze the potential in the 
Ekpyrotic/Cyclic model. Using numerical techniques we study the evolution of the moduli 
field focusing in particular on the transition from the contracting branch to the expanding 
branch. We show that a non-singular turn around of the scale factor along with the same 
for the moduli field is possible in the effective 4-dimensional description without any extra 
inputs. We summarize the results and discuss open issues in the concluding section. 



Such a modification to the Cyclic model potential has also been suggested in Refs. 361 137|. 
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II. BIANCHI-I MODEL: CLASSICAL THEORY IN ASHTEKAR VARIABLES 



The Bianchi-I spacetime is one of the simplest examples of spacetimes with anisotropies. 
It has vanishing intrinsic curvature and unlike Bianchi-IX model, the classical dynamics 
does not exhibit BKL behavior as the singularity is approached. The isotropic limit of this 
spacetime is k = FRW cosmological spacetime. We consider a homogeneous Bianchi-I 
anisotropic spacetime with a manifold X x R where £ is topologically flat (M 3 ). The spatial 
manifold is non-compact and in order to define the symplectic structure and formulate a 
Hamiltonian theory, it is necessary to introduce a fiducial cell V. The cell V has fiducial 
volume V = l\l 2 l 3 with respect to the fiducial metric q a b endowed on the spatial manifold. 
Here U refer to the coordinate lengths of the each side of the fiducial cell. 3 

Due to the homogeneity of the Bianchi-I spacetime, the Ashtekar variables take a sim- 
ple form. The matrix valued connection A l a and triad Ef can be expressed as Cj and Pi 
respectively (where i = 1,2,3) [38j|. The canonical conjugate phase space variables satisfy 

{ci,pj} = 8irGj8ij (1) 

where 7 = 0.2375 is the Barbero-Immirzi parameter. The triad pi are related to the three 
scale factors aj of the Bianchi-I metric 

ds 2 = —N 2 dt 2 + a\ dx 2 + a\dy 2 + a\dz 2 (2) 

as 

\p x \ = l 2 l 3 a 2 a 3 , \p 2 \ = /i/ 3 aia 3 , \p 3 \ = l 2 l 3 a 2 a 3 (3) 

where the modulus sign arises because of orientation of the triads (and is suppressed in the 
following) . 

The only non-trivial constraint to be solved in this model is the Hamiltonian constraint 
which when expressed in terms of Cj and pi takes the form: 

iV 

^cl = - „ 9T , (ciPi c 2 p 2 + c 3 p 3 c x p x + c 2 p 2 c 3 p 3 ) + T-C mSLt t (4) 

where H ma tt is the matter Hamiltonian which may describe perfect fluid and/or scalar fields. 
Physical solutions are obtained by the vanishing of the Hamiltonian constraint: 7i c \ ~ 
0. Equations of motions for the phase space variables are obtained by solving Hamilton's 
equations: 

pi = {pi, H c \} = -8nG-f — (5) 

and 

ti = {c h H c \} = 8nG-f — . (6) 
If we choose the lapse function N = 1, then using (0J and (J6]) we obtain 

d = 'ykdi = 'yk Hidi (7) 



3 As should be expected from any consistent treatment, physical predictions of this model are insensitive to 
the choice of k and one could choose these to be unity. We do not restrict to this choice in order to stress 
the independence of physics in the classical and especially the effective theory on the fiducial cell. The 
latter feature is not shared by alternative quantizations of Bianchi-I model in loop quantum cosmology. 
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where Hi = di/oj is the Hubble rate in the i-th direction. 

Similarly we can obtain the dynamical equations for the matter degrees of freedom. As 
an example, in case 7Y ma tt corresponds to a minimally coupled scalar field <fi with momentum 

(satisfying {</>,£></,} = 1) and potential V(4>), the dynamical equations take the standard 
form 

. d .3 

T^matt and p^ — — — W ma tt • (8) 



Taking the second derivative of <fi and using p^, we obtain the standard Klein-Gordon equa- 
tion for 0: 

+ ^^0 = -^V(0). (9) 

i 

Before we study the properties of dynamical equations it is useful to note the behavior of 
the classical Hamiltonian constraint under one of the underlying freedoms in the framework. 
This is to do with the change of the shape of the fiducial cell. Let us consider this change 
as: (h,h,h) — > (Ti?^^)- Under this change V — > V-J! 2 V 3 V and c 1; pi (and similarly other 
components) transform as 

Ci -> c[ = l[ci , pt -> p\ = \%p x . (10) 

Thus, the gravitational part of the constraint transforms as H grav — > l[l' 2 l' 3 'H gT a, v . It can 
be shown that the matter part of the constraint also transforms in the same way. Thus, 
^ciass/^ is invariant under arbitrary change in shape of the fiducial cell. 

Solving the classical constraint TC C \ ~ 0, dividing by the total volume V = V (aia 2 a 3 ) = 
(pip 2 Ps) 1//2 and using ([7]) we obtain the following equation relating directional Hubble rates 
with energy density 

H X H 2 + H 2 H 3 + # 3 #i = 8nG = 8nG p (11) 

where p is the energy density of the matter component: p = H mau tt/V. As expected, using 
( TTUj) one finds that Hubble rates and energy density are invariant under (li, l 2 , h) — > 1' 2 , 1's)- 

An interesting property of equations of motion for pi and Cj arises for matter with van- 
ishing anisotropic stress. It can then be shown that (pjCj — PjCj) is a constant of motion 
satisfying j38|, i.e. 

PiCi - pjCj = V(Hi - Hj) = ^V aij (12) 

where is a constant antisymmetric matrix. 

The directional Hubble rates can be considered to be the diagonal elements of an expan- 
sion matrix H^. The trace of this matrix is related to the expansion rate of geodesies in 
this spacetime as 

61=1(^+^ + ^3) = — (13) 
3 a at 

where a identified as the mean scale factor for our anisotropic model: 

a = (a^a^ 1 ^ . (14) 
The trace-free part of the expansion matrix leads to the shear term a if 

Oij = Hij - 05ij (15) 
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and defines the shear scalar a 2 = o^a^ given by 

3 1 
J2(^- ) 2 = o ((Hi ~ H 2 f + (H 2 - H 3 f + (H 3 - H,) 2 ) 



a 2 



3 

= ("12 + "23 + «3l) ( 16 ) 

where to obtain the last expression we have used Eq. (TT2]) . 

The generalized Friedman equation for the mean scale factor can be obtained by consid- 
ering the mean Hubble rate 

a 
a 



H 



^X>, (17) 



which yields 

H 2 = ^(H 1 H 2 + H 2 H 3 + H 3 H 1 ) + -L _ H 2 f + (H 2 - H 3 ) 2 + (H 3 - H^ 2 ) . 
On using Eqs.lfTTT) and ( fTBT) we obtain 

H= — " + ^- < 18) 

with 

S 2 = iaV. (19) 

From f|T6|) it follows that the shear scalar S is a constant of motion in the classical theory. 
As we will show later, this feature does not hold in the loop quantization where only at the 
classical scales E approaches a constant value. 

Analysis of the generalized Friedman equation (fTSj) immediately shows that the singularity 
is inevitable in the classical theory as the scale factor approaches zero. It is important to 
note that in the classical theory anisotropic term dominates on approach to singularity 
unless p corresponds to matter with equation of state w > 1 i.e. with an equation of state 
of an ultra-stiff fluid. In the case when matter has equation of state w < 1, the energy 
density grows slower than anisotropic term as a(t) — > 0. We will discuss later that in the 
Ekpyrotic/Cyclic model the equation of state w ^> 1 during the ekpyrosis phase (when the 
moduli is in the steep negative region of the potential). This causes anisotropic term to be 
subdued in subsequent evolution. Thus, as the universe approaches classical big bang/crunch 



singularities in this model one expects anisotropies to become very small [39 



III. EFFECTIVE DYNAMICS 

In the loop quantization the elementary variables are the holonomies of the connection 
and the flux of the triad (related by a constant for the present case). Elements of the 
holonomies are of the form exp(z/iCj) where fi labels the edge length along which a holonomy 
is computed. The classical Hamiltonian constraint is expressed in terms of holonomies and 
triads and then quantized. The resulting quantum constraint is generically non-singular and 
the picture at the Planck scale in Bianchi-I model turns out to be similar to the one in the 
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isotropic model |24J. One can also derive the effective Hamiltonian using similar techniques 
as used for the isotropic model [3l[ and for iV = 1 is given by 4 

If sin(^iCa) sm(/z 2 c 2 ) \ 
ries = - „ 9T , = z V\V2 + cyclic terms + 7Y matt (20) 

where 



/ii = AW— , /! 2 = A A /-^-, and p 3 = A,/— . (21) 
V V PiKs V 

Here A arises due to regularization of the field strength of the connection by the underlying 
quantum geometry. Its value is given by 

A 2 = 4V3tt74j . (22) 

It is to be noted that the effective dynamics resulting from the above Hamiltonian is invariant 
under the choice of the fiducial cell. In particular this is the only known loop quantization 
of Bianchi-I spacetime which is independent of the shape of the fiducial cell [30j. Comparing 
([20]) and (|3]) we notice that the change from classical to effective theory only comes in the 
form of replacement of connection components Cj with sin(/2/C/)//i/. Under the transforma- 
tion: (li, l 2 , 13) — > 1' 2 , 1' 3 ), fii — > fl'j = (l/l'j)jjLi. Hence sm(fijCj) / fij transforms in the same 
way as C/. Thus transformation properties remain same as in the classical theory and the 
resultant dynamics and physical predictions are unaffected by the freedom of the choice of 
the fiducial cell. 5 

Substituting (1211) in (120]) and solving for the Hamiltonian constraint 7i e ff ~ we obtain 
P = o ^2x2 ( sin G"i c i) sin(/i 2 C2) + cyclic terms) . (23) 

o7T(jr7 A 

Since sin(/2jCj) are bounded functions, above equation implies that the energy density can 
never diverge in the effective loop quantum cosmology. The maximum value of the terms in 
the parenthesis determines the upper bound for the energy density: 

3 

P — Pmax; Pmax = g^Q^y2^2 ' \ ) 

The upper bound for the energy density in the Bianchi-I anisotropic model coincides with 
the value in the isotropic model 0]. 

It should be noted that the matter Hamiltonian 7i mat t in fl20l is unmodified from its 
classical expression, due to which the dynamical equations for matter remain the same 



4 In order to compare with earlier works (for eg. equations in Ref.[30l.l38j]). in our convention \ii correspond 
to often used fi^. 

This is in contrast to other proposals for quantization of Bianchi-I spacetimes such as in Ref.kjJ and 



those motivated by lattice refinement considerations [42J. As emphasized in Ref. [30(, since the effective 
dynamics in these models is not invariant under the change in shape of the fiducial cell, they lack predictive 
power. It turns out that in the effective dynamics of these proposals the expansion factor and shear are 
also unbounded, even if one fixes the fiducial cell. 
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as Eqs.® and (Q. Modified dynamical equations can be obtained from the Hamilton's 
equations using (|2Q|) . As an example, for the triad component p\ we get 



Pi 



-8ttG 
Pi 



dci 



— cos(;UiCi) (sin(/i 2 c 2 ) + sin(/i 2 c 2 )) 
7A 



(25) 



(and similarly for p 2 and p 3 ). 

In order to find the equation for directional Hubble rates we first note that from Q we 

get 



h 



P2P3 
Pi 



1/2 



a-2 



1 



P-iPl 
P2 



1/2 



and 03 



1 

h 



P1P2 

P3 



1/2 



(26) 



Taking their derivatives and using Eq.( [25l) and corresponding equations for p 2 and ^3 we 
obtain 

Hi = — = —— (sin(/iici - /2 2 c 2 ) + sin( / uici - ^3) + sin(/i 2 c 2 + fi 3 c 3 )) . (27) 
a\ 27A 

Similar equations can be derived for the directional Hubble rates if 2 and i?3. It is clear from 
the above equation, unlike in the classical theory the directional Hubble rates in effective 
loop quantum dynamics are always bounded. 

The Hamilton's equation for connection components can be derived in a similar way As 
an example, for c\ one obtains 



C\ = 87rG7 



dH, 



eff 



<9ci 



2pi7A 



c 2 p 2 cos(/i2C 2 )(sin(/iiCi) + sin(^ 3 c 3 )) + c 3 p 3 cos(/x 3 c 3 )(sin(/xiCi) + sin(/i 2 c 2 )) 



cipx cos(/i 1 ci)(sin(/i 2 c 2 ) + sin(// 3 c 3 )) - H1P2P3 



sin(/x 2 c 2 ) sin(/i 3 c 3 ) 



+ sin(/iiCi) sin( / u 2 c 2 ) + sin(/i 3 c 3 ) sin^ci) 



+ 8vrG7 



P2P3 / P 



Pi 



dp 
2 dpi 



(2f 



Using (1251) with ( 1281) and similar equations for othe connection and triad components, it 
can be shown that matter with vanishing stress satisfies 



d 
It 



{PiCi - PjCj) = 



(29) 



Thus as in the classical theory (jp constant of motion in the effective loop 

quantum dynamics. However, due to Eq. (l2"7j) . unlike in the classical theory the relation (J7J) 
is no longer satisfied. A consequence is that in the effective dynamics 



p^ - pjCj ^ V{Hi - Hj) . 



(30) 
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In the classical approximation //jCj <C 1, Eg. (1271) leads to q ~ 7/idj which implies 
PiCi — pjCj ~ V(Hi — Hj) as we approach the classical scales. 

We can now evaluate the shear scalar a 2 in the effective description of LQC Since (Hi — 
Hj) are no longer constant except at the classical scales, we note from Eq. (ll6l) that a 2 is not 
a constant in the effective LQC. Using the Hamilton's equations for p^, a straightforward 
calculation leads to 



a 2 



1 



cos(^ 3 c 3 )(sin( / u 1 ci) + sin(/i 2 c 2 )) - cos(/i 2 c 2 )(sin(/iiCi) + sin(/i 3 c 3 )))' 
+ (cos( / u 3 c 3 )(sin( / u 1 ci) + sin(jU 2 c 2 )) - cos(/2iCi)(sin(p 2 c 2 ) + sin(/i 3 c 3 ))) 2 
+ (cos(/2 2 c 2 )(sin(/iiCi) + sin(/2 3 c 3 )) - cos(/ZiCi)(sin(/2 2 c 2 ) + sin(/i 3 c 3 ))) z 



(31) 



This implies that S 2 defined via Eq. (TT9T) is not a constant of motion in the effective loop 
quantum dynamics. Nevertheless, it is clear that a 2 is bounded above by a fundamental 
value in loop quantum cosmology. A detailed analysis of the behavior of anisotropy and 
energy density will be done else where. It is worth pointing out some of the notable features 
of the effective dynamics. These include: 

(i) Unlike the bounce in the isotropic model, due to presence of anisotropies the energy 
density at the bounce may be less than p max . In fact it is possible for the bounce to occur 
with p = 0. In this case bounce occurs purely because of interplay of quantum geometric 
effects with anisotropy. 

(ii) The maximum value of anisotropic shear at the bounce is 

^max = ■ (32) 

(iii) It is also possible for the effective dynamics to saturate Eq. (1241) at the bounce. However, 
in that case anisotropy vanishes at the bounce. 

These features exhibit richness of the effective dynamics in the Bianchi-I model. They 
result due to incorporation of quantum geometry effects in the effective Bianchi-I spacetime. 
When the components of the spacetime curvature become large, there are significant depar- 
tures from the classical theory. When these components are small compared to the Planck 
scale then dynamical equations approximate their classical counterparts and one recovers 
the classical description. It is important to note that the upper bounds on energy density, 
shear and Hubble rate are direct consequences of the underlying quantum geometry. In the 
classical limit A — > 0, p ma x, -f^max and <r max diverge and the evolution is singular. The bound- 
edness of these quantities plays the crucial role to construct non-singular Ekpyrotic/Cyclic 
model as we we discuss in the following section. 



IV. EVOLUTION WITH THE CYCLIC POTENTIAL 



In the previous section we analyzed the effective dynamics of Bianchi-I model in loop 
quantum cosmology for matter with vanishing anisotropic stress. An important feature of 
quantum geometry modified dynamics turns out to be singularity resolution. We showed that 



10 



the energy density, Hubble rates and shear are always bounded in loop quantum cosmology 
independent of the equation of state of the matter content. This result is very encouraging 
for the primary question posed in this work: Is there a viable non-singular evolution with 
the Ekpryotic/Cyclic potential in the effective loop quantum cosmology? We now explore 
the answer to this question. 

Let us first briefly recall some of the salient features of Ekpyrotic/Cyclic model which is 
motivated by M-theory. It is hypothesized that the observable universe is constrained on 
a visible brane which interacts gravitationally with a shadow brane. The collision between 
these branes constitutes the big bang/big crunch singularities for the observable universe 
in the 4-dimensional effective description. The inter-brane separation is determined by a 
moduli field ((f)) and the interaction potential is given by 



V = V (l - e^) exp(-e^) (33) 

where V a , o\ and 02 are the parameters of the potential. 

An attractive feature of the Ekpryotic/Cyclic models is the different roles which the 
moduli plays in various epochs. When the branes approach each other, the moduli field 
moves slowly in an almost flat and positive part of the potential which in the effective 4- 
dimensional description leads to a period of dark energy domination. Eventually when the 
potential energy of the moduli balances the Hubble expansion rate, the moduli field turns 
around and the inter-brane separation starts decreasing leading to a contracting phase in the 
universe. When the branes approach each other, the moduli potential is steep and negative 
leading to an ultra-stiff equation of state, i.e. w > 1. It turns out that this feature leads 
to a decrease in the anisotropy of the universe at small scale factors [39(. The ultra-stiff 
equation of state is also responsible for producing a scale invariant spectrum of fluctuations 
with a significant non-gaussian contribution (43[. The latter serves as a distinct prediction 
to the inflationary scenarios, and is available for test in the near future. The branes collide 
at = —00 which follows their separation and running of moduli field towards the positive 
value of potential, leading to a cyclic phenomena. Thus, the 4-dimensional effective dynamics 
description of the Ekpyrotic/Cyclic model is very rich and interesting. 

A primary issue in the Ekpyrotic/Cyclic model is to understand the transition from the 
contracting to the expanding branch. This transition must be non-singular and is vital to 
propagate the perturbations generated prior to the collision of the branes to the regime 
following the collision. In the 5 dimensional picture, the bulk spacetime near the collision of 
the branes appears as a compactified Milne space and it has been argued that the resulting 
singularity is milder than the conventional big bang/ crunch singularity [44j|. Further, in this 
picture it may be possible to match the perturbations across the transition [13J . This however 
results in the mixing of the modes with a strong sensitivity of predictions to the process 
of transition which ignores quantum gravitational effects. Other ideas to understand this 



transition have been proposed including perturbative string theoretic treatment [45| and 



AdS/CFT correspondence [14j, Il5l |. Though promising directions to explore, these reveal 
little about the transition from the contracting to the expanding branch in the 4-dimensional 
effective theory. If one assumes that the 4-dimensional effective theory has no inputs from 
a theory beyond GR, this difficulty is not surprising. The reason is that in GR such a non- 
singular turn around is forbidden unless matter violates null energy condition. Though it 
is not possible to have this violation in the Ekpyrotic/Cyclic models, a variant of these has 



been proposed in form of New Ekpyrotic model [46|, |47| which includes a ghost condensate. 



The latter is responsible for violation of null energy condition and can result in a bounce. 
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Whether such a non-singular bounce is a generic feature of this model and is not affected 
by the instabilities is not an open issue. Further, it is not clear whether the 4- dimensional 
effective picture with a ghost condensate as conjectured in this model is valid at high energies 



of interest 48 



Given that non-perturbative quantum geometric modifications and their role in singular- 
ity resolution is well understood in the homogeneous models of loop quantum cosmology, 
we now explore the effective 4-dimensional loop quantum dynamics of the moduli field to 
understand the transition from the contracting to the expanding branch. An underlying 
assumption of such an analysis is the treatment of potential fl33l) as an effective one for the 
field in four spacetime dimensions. This enables us to use effective equations derived in 
the previous section with as the source of matter. An investigation on these lines was 



carried out earlier with an additional assumption that the spacetime be isotropic (35|. Its 
main conclusions were: 

(i) the big bang/big crunch singularity is generically avoided irrespective of the choice of the 
parameters of the potential and initial conditions for the field 0. 

(ii) However, it is not possible in the isotropic effective loop quantum description for the 
field to turn around from = — oo. Thus a cyclic evolution as envisioned in the Ekpy- 
rotic/Cyclic model paradigm was not possible in the isotropic loop quantum cosmology. 6 

It should be noted that the absence of the turn-around of from large negative values 
is not a limitation of the loop quantum dynamics, but is a feature shared also by the 
classical theory. This can be easily seen from the classical Friedman equation for the isotropic 
model. For simplicity let us assume that the only contribution to the matter energy density 
originates from the moduli field 

H 2 = — P<i> (34) 



where 



~ 2 + V{4>) . (35) 



For the moduli field to turn around it is essential that must vanish. Using ( 135]) in ( 1341) we 

get 

It is clear that a turn around of the field is not possible when the potential is negative. 
Since in the Ekpyrotic/Cyclic model the potential is negative (and asymptotes to zero) 
when the branes are supposed to collide, can not vanish and the moduli can not turn 
around. Presence of additional energy density components such as radiation is unlikely to 
affect this conclusion. The reason being that in the vicinity of brane collision (a — > 0), 
the dominant contribution to the energy density comes from the kinetic energy of the moduli. 

Interestingly, above situation changes if the spacetime has a non-vanishing anisotropic 
shear. In this case, the classical generalized Friedman equation yields 



6 If the inter-brane potential is modified such that it is positive for large negative values of <j>, then the 



turn around of the field and a cyclic description is possible 35] . A modification to the cyclic potential as 
proposed in the bicyclic scenario [3^ would hence result in a desired non-singular evolution. 
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(37) 



The shear term can be seen as an additional component of energy density which scales the 
same way as p^. Thus depending on the strength of the shear term, it is possible that <fi 
can turn around even when V{4>) < 0. However, the classical theory is singular and the 
above equation breaks down when spacetime curvature becomes very large. Hence, though 
we get some insights on the possible role of anisotropic to alleviate the above problem, it 
is not possible to obtain a non-singular cyclic evolution without any additional inputs from 
quantum gravity or going beyond the scope of the classical theory. 

Given the non-trivial role the shear term can play in the dynamics of the moduli field 
in this model, it becomes important to analyze the effective loop quantum equations with 
anisotropies. For simplicity we assume that the only contribution to the matter density 
originates from the field 0. By specifying the 7Y ma tt corresponding to the Cyclic potential 
jil we can derive the dynamics from Hamilton's equations for q, pi, the moduli </> and its 
momentum p^ as obtained in Sec. III. These equations can be numerically integrated to 
obtain the behavior of directional scale factors, the mean scale factor, Hubble rates, energy 
density and shear (which is calculated from its definition (j!6p ). 




FIG. 1: Plot of the mean scale factor (in Planck units), a = (aia2as) 1 ^ is shown for the initial 
conditions (all values in Planck units) <fi = 0.43, </> = —0.038, p± = 64, p2 = 72, p% = 68, c% = —0.8, 
C2 = —0.7 and S 2 = 5.80407. The mean scale factor experiences multiple small bounces as a result 
of bounces of the individual scale factors (aj). Unlike the classical theory, there is a non-singular 
transition from the contracting to the expanding branch. 

We performed various numerical simulations with different initial data for the moduli 
field, expansion rates and anisotropy. The initial conditions to solve the dynamics were 
provided for a contracting universe with small initial anisotropy and the value of the field <fr 
in the positive part of the Cyclic model potential. Since the primary aim of our analysis is 
to investigate the resolution of singularity, we choose without the loss of generality param- 
eters and initial conditions independent of the consideration from the ones constrained by 
observations. We now illustrate some numerical results obtained from the effective dynamics. 
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FIG. 2: Evolution of the moduli field in Planck units is shown for the initial conditions in Fig. 1. 
The moduli field starts from a positive region of the potential, rolls to the negative part and turns 
around in the Planck regime. After the turn around the moduli again reaches the positive part of 
the potential, ready for another cycle. 
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FIG. 3: The plot shows the variation of shear S 2 in the transition region from contracting to 
expanding branch for initial conditions in Fig. 1. Considerable variation for a short period before 
and after the bounce is evident from the spikes in value of log(S 2 ). We also see that in the classical 
regimes at both small and large values of t, shear approaches similar values. 



The first set of evolution is depicted in Figs. 1-4. Parameters in the potential fl33|) 
were chosen as V a = 0.02, o\ = 0.3v^8~7r and a 2 = COOv^S^- (For simplicity we choose 
G = h = c = 1). Initial conditions for the moduli field were = 0.43 and <fi = —0.038. Initial 
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FIG. 4: Variation in directional Hubble rates is plotted for the initial data in Fig. 1. Green curve 
corresponds to Hi, red to H2 and the blue curve corresponds to H%. The spikes reflect bounces 
and recollapses of individual scale factors in the Planck regime during the transition from the 
contracting to the expanding branch. 

conditions (provided at time t = 0) for the triad and connection components were pi = 64, 
p 2 = 72, p 3 = 68, c\ = —0.8, c 2 = —0.7 with the initial anisotropy £ 2 = 5.80407. (The 
connection component C3 was determined by solving the Hamiltonian constraint 7i e K ~ 0). 
These initial conditions correspond to a contraction of all the scale factors, i.e. di < 0, a? < 
and 0,3 < 0. Fig. 1 shows the evolution of mean scale factor a = (a^a^) 1 ^ 3 . We find that the 
mean scale factor undergoes multiple bounces and recollapse in the Planck regime and there 
is a non-singular evolution from the contracting to the expanding branch. This in a sharp 
contrast to the classical theory where such a transition is forbidden and the contracting and 
expanding branches are disjoint in the 4-dimensional description. Interestingly as depicted 
in the Fig. 2 we also obtain a turn around of the moduli field. The moduli starts from the 
positive part of the potential when the universe is contracting, rolls down the negative region 
of potential and when the spacetime curvature reaches Planck value, it turns around and 
goes back to the positive part. In the subsequent evolution the classical friction term stops 
the moduli field and it turns around causing a contraction of the universe. Thus leading 
to a cyclic phenomena. In the above numerical simulation, the turn around of the moduli 
field occurs at approximately t ~ 100 (in Planck units) coinciding with the midpoint of the 
transition from contraction to expansion of the scale factor. 

Fig. 3 shows the evolution of the shear term obtained by the numerical integration. We 
find that the shear term varies significantly from its classical value in the region where the 
transition occurs. However, it is everywhere bounded and does not affect the occurrence of 
bounces. As discussed earlier, due to the property of evolution of shear in loop quantum 
cosmology, the value of shear scalar in the low curvature regime before and after the 
transition turns out to be the same. Also for this numerical run, the shear scalar is very 
small when the moduli field turns around. We also find that the behavior of S 2 consists of 
spikes in the Planck regime where the mean scale factor of the universe undergoes multiple 
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bounces. Using the expression of shear factor (obtained from (fT9j) and (fl6l) for loop quantum 
evolution), the cause of these spikes turns out to be due to rapid variation in the Hubble 
rates in different directions occurring during multiple bounces. This is illustrated in Fig. 4 
where in different curves we have shown the behavior of directional Hubble rates. As can be 
seen the spikes in shear scalar coincide with the spikes in the directional Hubble rates which 
occur due to multiple bounces and recollapse of anisotropic scale factors in the Planck regime. 




FIG. 5: Plot shows evolution of the scale factor in the presence (solid curve) and absence (dashed 
curve) of anisotropy for initial conditions 4> = 0.4, <f> = —0.03, p\ = 64, p 2 = 72, p 3 = 68, c\ = —0.6, 
C2 = —0.5 and £ 2 = 9.2365. The classical singularity is avoided in both the cases. 



Another example of results obtained from numerical integration are shown in Figs. 5-8. 
The parameters of the potential (I3"3"j) were chosen as V a = 0.01, a% = 0.3y/8n and a 2 = 
0.09v87r. The initial conditions specified at time t = were <p — 0.4, = —0.03, p\ = 64, 
p 2 = 72, p 3 = 68, C\ = —0.6, c 2 = —0.5 and £ 2 = 9.2365. These initial conditions correspond 
to two of the anisotropic directions contracting (dj < 0, 0,2 < 0) and one expanding (0,3 > 0). 
To understand the non-trivial role played by anisotropic term we also performed analysis 
of evolution in the isotropic loop quantum model for the same initial conditions for moduli 
field and the isotropized initial conditions for triad and connection. In Fig. 5 we have shown 
the evolution of the scale factor in the present model (solid curve) and the isotropic model 
(dashed curve). We find that the bounce of the scale factor is present both in presence 
and absence of anisotropy. The loop quantum evolution with the Cyclic potential is thus 
non-singular with or without anisotropy. Fig. 6 shows the evolution of the moduli field and 
we find that as for the previous case (Fig. 2) there is a turn around of the moduli field from 
the negative values of the potential to the positive part resulting in a cyclic phenomena. As 
before, the turn around of the field occurs at the midpoint of the transition of the mean 
scale factor from the contracting to expanding branch. Fig. 6 also shows that for the same 
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FIG. 6: Evolution of the moduli field is compared for initial conditions in Fig. 5 for the anisotropic 
(blue) and isotropic (red) model. As in Fig. 2 the moduli field turns around in the presence of 
anisotropies in the Planck regime and leads to a non-singular cyclic model. It fails to turn around, 
and approaches (f> = — oo, in the absence of anisotropies in confirmation with results of Ref. 35]. 
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FIG. 7: The shear scalar E is plotted for the numerical run in Fig. 5. We see that its value is 
preserved before and after the transition period. The shear term is bounded across the evolution 
with spikes corresponding to multiple bounces and recollapses of the individual scale factors. 



initial conditions the turn around of the moduli is absent for the isotropic model. Thus, 
anisotropies play a very important role to obtain a non-singular viable cyclic description. 

The evolution of shear and energy density are depicted in Figs. 7 and 8. We see that 
the shear term undergoes a significant variation in the region of non-singular transition but 
remains bounded. The spikes in the variation correspond to the spikes in the directional 
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FIG. 8: Plot of the energy density of the moduli field in the transition regime is shown. Bounces 
of the individual scale factors result in spikes in its behavior. Due to presence of non-vanishing 
anisotropy, the energy density at these bounces does not need to saturate to its maximum allowed 
value. 




FIG. 9: Plots of the mean scale factor a and the scalar field <j) are shown for the initial conditions 
with very small anisotropics. These are <f> = 0.4, <ft = —0.03, p\ = 64, p2 = 72, p3 = 68, c\ = —0.6, 
c 2 = -0.532 and S 2 = 0.8517. 

Hubble rates which occur due to multiple bounces and recollapses of the scale factors 
(01,02,03). Further, its value before and after the transition is preserved. Unlike the 
previous case in Fig. 2, the shear scalar is not small at the turn around of the moduli field 
but shows a spike. The behavior of energy density is shown in Fig. 7 which demonstrates 
that at the bounces of the directional scale factors (associated with the occurrence of 
spikes) the energy density does not saturate to its maximum value, a feature due to the 
non- vanishing anisotropy at bounces. In contrast, the energy density at the bounce in the 
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FIG. 10: This plot shows the evolution of different scale factors for the simulation corresponding 
to Fig. 9. We have plotted the logarithm of directional scale factors. The green curve (hidden 
behind red curve due to small anisotropy between a\ and a>2 directions) depicts a\, red curve a 2 
and blue curve 03. 




FIG. 11: These plots depict evolution of scale factor and scalar field for initial conditions: = 0.4, 
j> = -0.03, pi = 64, p 2 = 72, p 3 = 68, cy = -0.6, c 2 = -0.1 and S 2 = 1449.4278. Unlike the 
previous cases, we find that the asymmetry in the mean scale factor and the moduli field around 
the bounce point is enhanced. Note that the moduli field reaches the positive value for t ~ 220 
before starting a new cycle. 



isotropic model always reaches its maximum value, p cr it = 0.41ppi at the bounce. 

In both of the previous simulations, though the initial anisotropic shear is small it is 
nevertheless of the same order. An interesting question which arises is whether the non- 
singular transition of the scale factor and the turn around of the moduli field occurs if 



the shear is decreased. It is to be noted that as shown in Ref.[30| the loop quantization of 



Bianchi-I model leads to an upper bound on the expansion rate and anisotropic shear. These 
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bounds are generic and independent of initial conditions for the matter content. The upper 
bound on the expansion rate implies that the mean scale factor would always bounce for the 
Cyclic model potential as the classical singularity is approached. Further, from the discussion 
of the classical theory we recall that the turn around of moduli field would require appropriate 
choice of initial conditions of the scalar field and anisotropy. In various simulations which 
we performed, we found that there exists a large range of initial data for which the turn 
around of moduli field can occur even for very small anisotropies. Results from one such 
example are depicted in Figs. 9-10. We chose the same values of parameters as in previous 
simulations. The initial values are p\ = 64, p 2 = 72, p 3 = 68, ci = — 0.6, c 2 = —0.532 and 
X 2 = 0.8517. As can be seen from these plots a significant decrease in the initial anisotropy 
does not affect the bounce of the scale factor or the turn around of the moduli field. 

It is interesting to note that in these simulations the behavior of the mean scale factor 
and the scalar field seems symmetric across the middle of the transition from the contracting 
to the expanding branch. To understand this we first note that the weak symmetry of 
the variation of the scalar field in the transition regime stems from lack of interaction 
with any other form of matter. We expect that presence of additional matter degrees 
of freedom would lead to a pronounced asymmetry. Further note that in the plots we 
depict the mean scale factor which actually suppresses the asymmetry of the individual 
scale factors across the bounce. This becomes clear if we plot the anisotropic scale factors 
as shown in Fig. 11 for the numerical simulation for Fig. 10. In Fig. 12 we show 
results from another numerical simulation where the asymmetry in the mean scale factor 
is not suppressed. The parameters are chosen as same in previous cases and the initial 
conditions are p 1 = 64, p 2 = 72, p 3 = 68, c : = -0.6, c 2 = -0.1 and £ 2 = 1449.4278. As 
we can see the behavior of scale factor and the scalar field across the bounce is not symmetric. 

We summarize the main results of the numerical analysis as follows: 

1. Irrespective of the choice of initial conditions, the classical singularity at a = is 
generically avoided in effective dynamics of loop quantum cosmology for the Cyclic 
model potential. 

2. Starting from arbitrary anisotropic conditions, scale factors in different directions 
bounce when spacetime curvature becomes close to the Planck value. This causes 
a bounce of the mean scale factor in the Planck regime. The shear term remains 
bounded during the evolution and it approaches the constant classical value when 
spacetime curvature becomes small. 

3. Interestingly, it is not difficult to choose the initial conditions such that avoidance 
of the singularity is accompanied by a turn around of the moduli field in the 
negative regime of the potential. The moduli then rushes towards the positive 
part of the potential, stops at certain positive value of the potential and roll back 
towards the negative value of the potential. This leads to a cyclic model of the universe. 
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V. SUMMARY AND OPEN ISSUES 



Ekpyrotic/ Cyclic model is a very interesting paradigm for the early universe which is 
considered as an alternative to the inflationary scenarios. A key issue in this model is to 
understand the transition from contracting to the expanding branch in the 4-dimensional 
spacetime picture. In the 5-dimensional picture this transition corresponds to the collision 
between boundary branes in the bulk. Though novel insights have been gained in the latter 



phenomenon [14J, [15|, |44J, |45| , obtaining a non-singular transition in the 4-dimensional picture 
has remained an open issue. Understanding which is important for various reasons including 
the way cosmological perturbations propagate from the contracting to the expanding branch. 
Given the Penrose-Hawking singularity theorems, a non-singular transition is not possible 
in the framework of classical theory. There are hopes to alleviate this problem in the 
classical framework by introduction of a ghost condensate, however the approach has its 



own limitations 46, 48 



Our analysis is based on the widely accepted notion that resolution of singularities involves 
going beyond the classical description of gravity. In particular, inputs from non-perturbative 
quantization of gravity may be necessary to avoid singularities. In our approach we have used 
the effective spacetime description of loop quantum cosmology to analyze the dynamics of the 
Ekpyrotic/Cyclic model. It is a non-perturbative quantization of homogeneous spacetimes 
based on loop quantum gravity and has successfully addressed resolution of cosmological 
singularities in various settings j3, @, 12, H, 24] with a general picture of replacement of big 



bang with big bounce at the Planck scale. In this work we have focussed on the resolution of 
big bang/crunch singularity in the Ekpyrotic/Cyclic model. Such an investigation has been 
performed earlier [35( 7 using the assumption that the spacetime be purely isotropic. It was 
found that though the 4-dimensional scale factor generically bounces in the Planck regime 
a viable Cyclic model is not possible due to lack of a turn around of the moduli field from 
the negative region of the Cyclic model potential. 



To probe whether conclusions reached in the previous work [35( were artifacts of ignoring 
anisotropic shear, we include them in the analysis and investigate the dynamics with the 
effective description of loop quantization of the Bianchi-I spacetime. Our analysis assumes 
vanishing intrinsic curvature (whose consequences are discussed below) and the moduli field 
as the only source of the matter energy density. The effective dynamical equations are 
complicated to solve analytically due to which various numerical simulations were performed 
with different choices of the parameter of the Cyclic model potential and initial conditions. 
We find that the universe undergoes a non-singular transition from the contracting to the 
expanding branch accompanied by small multiple bounces for individual scale factors in 
the Planck regime. In contrast to the previous results, we find that it is possible to easily 
choose initial conditions such that the moduli field turns around from the negative part 
of the potential in the Planck regime. Thus leading to a potentially viable non-singular 
Ekpyrotic/Cyclic model without introduction of any exotic matter. This is the novel result 
of our analysis. 

It should be noted that the anisotropies play a very non-trivial role to obtain a non- 
singular cyclic model in our analysis. At first sight it seems perplexing because in the 
ekpryosis phase, which occurs when the moduli field is in the steep negative region of the 



7 See also Ref. [49| for an earlier work which ignored modifications to the gravitational part of the Hamil- 
tonian. 
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potential, anisotropies become small and the universe evolves towards an isotropic phase. 
However as we demonstrated, even though anisotropies may become very small during the 
evolution, their non-zero value is important to turn around the moduli field in the negative 
part of the potential. Recall that such a turn around is not possible even for the classical 
theory in the absence of anisotropies. Hence the loop quantum evolution successfully leads 
to a turn around of the scale factor and the moduli without affecting nice features of the 
Ekpyrotic/Cyclic model. 

There remain several open questions which require further investigations. One of them 
deals with generalizing our model to include the intrinsic curvature of the spacetime which 
will require studying loop quantization of spacetimes such as Bianchi-IX. It will also help in 
obtaining insights on the BKL behavior in loop quantum cosmology and whether it affects 
the conclusion reached in this work. A second issue is to gain more analytical control on 
the effective dynamics of the anisotropic loop quantum cosmology. Due to its complexity it 
is difficult to obtain an analog of generalized Friedman equations with shear term in loop 
quantum cosmology, however various insights have been obtained on general features of the 
effective dynamics. These include, showing that expansion parameter and shear term are 
universally bounded j30|. Further analytical studies in this direction on the Ekpyrotic/Cyclic 
model will be reported elsewhere. It will be also useful to understand the full loop quantum 
dynamics with the Cyclic potential and compare with the effective dynamics treatment as 
done presently. Also in order to have a more realistic Ekpyrotic/Cyclic model, it is impor- 
tant to include radiation, study its interaction with moduli field and its influence on the 
non-singular transition obtained in this work. Finally, it is an open problem to introduce 
inhomogeneities in our framework and understand the role of quantum gravitational effects 
on their propagation through the bounce. Recent results on study of Fock quantized inho- 



mogeneous modes on the loop quantum spacetime are encouraging in this respect 50]. It 
is quite possible that incorporation of these inhomogeneities may reveal subtle imprints of 
quantum gravity on the predictions of the Ekpyrotic/Cyclic model. 
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